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In many electron spin qubit systems coherent control is impaired by the fluctuating nuclear spin 
bath of the host material. Previous experiments have shown dynamic nuclear polarization with 
feedback to significantly prolong the inhomogeneous dephasing time Tf by narrowing the distribution 
of nuclear Overhauser field fluctuations. We present a model relating the achievable narrowing of the 
Overhauser field to both the pump rate and the noise magnitude and find reasonable agreement with 
experimental data. It shows that former experiments on gated GaAs quantum dots were limited by 
the pump rate of the pumping mechanism used. Here we propose an alternative feedback scheme 
using electron dipole spin resonance. Sequentially applying two ac electric fields with frequencies 
slightly detuned from the desired Larmor frequency results in a pump curve with a stable fixed 
point. Our model predicts that Tf values on the order of microseconds can be achieved. 


One limiting factor of spin qubits based on semicon¬ 
ductor quantum dots is the fluctuation of the effective 
magnetic field arising from nuclear spins. The inter¬ 
action of these nuclear spins in the host material with 
the electron spin in the quantum dot induces decoher¬ 
ence. Several techniques have been developed to reduce 
decoherenceii— Narrowing the distribution of the nuclear 
fluctuations was shown to prolong the inhomogeneously 
broadened dephasing time for gate defined GaAs quan¬ 
tum dots^ and for self-assembled InAs quantum dotsi^i 
To reduce fluctuations the nuclear field is regulated by 
dynamic nuclear polarization (DNP). DNP controls the 
Overhauser field by transferring the electron spin angular 
momentum to the nuclear spins via the hyperfine interac¬ 
tion. If the effectiveness of the polarization step depends 
on the current state of the nuclear spin bath, the qubit 
can be used as a closed feedback loop. 

In gate defined GaAs double quantum dots, a roughly 
tenfold enhancement of T 2 * has been achieved by Landau- 
Zener-sweep driven electron-nuclear spin flip-flops that 
were preceded by a free evolution of the qubit to create 
a feedback effect £ In self-assembled InAs quantum dots 
a DNP feedback effect was also observed as the quan¬ 
tum dot resonance exhibited locking to frequencj^i^ or 
repetition rate^ of the incident laser. Driving GaAs gate 
defined quantum dots with electron dipole spin resonance 
(EDSR) showed locking to the spin-resonance condition 
as well— 

The narrowing effect has been examined theoretically 
in various studies but these were focussed on microscopic 
models and optical manipulation— - — A phenomenolog¬ 
ical rate-equation model to describe the effectiveness of 
DNP in reducing the Overhauser field variance was de¬ 
rived by Vink et a/— and by Latta et al~ The model 
relates the maximum nuclear field achievable by polar¬ 
ization to the nuclear field variance by treating the po¬ 
larization as a random variable influenced by DNP and 
relaxation proportional to the current value. This model 
was further refined by Yang et aZ— However, the domi¬ 


nant relaxation process in gated dots is spin diffusion of 
the local polarization to its surroundings, which is not a 
simple relaxation process : 17 : 18 As a result, the maximum 
field is determined by the slow long range nuclear dynam¬ 
ics, while the variance of the fluctuations also depends 
on short range fluctuations. Overall, diffusive dynamics 
have received little attention, but are of particular inter¬ 
est for gated dots. 

In this paper we present a model relating the achiev¬ 
able narrowing of the nuclear field distribution to the 
pump rate and the noise magnitude by explicitly consid¬ 
ering the different relaxation rates for polarization modes 
with different wave length. To this end we solve the spin 
diffusion equation including Langevin dynamics. The re¬ 
sults are in reasonable agreement with experimental data 
and show that the variance of nuclear fluctuations de¬ 
pends directly on the pump rate via the feedback gain. 
Previous experiments were limited by the pump rate, for 
example in gate defined GaAs dots the pump rate was 
set by slow Landau-Zener sweeps at the avoided cross¬ 
ing of the singlet and T + state (’ST+ scheme’)?^ Thus, to 
overcome the limitation of the ST + scheme, we propose a 
different polarization scheme based on the work of Laird 
et aZjiS the ’EDSR scheme’. Here an ac electric field is 
used to flip the spin of the electron in the dot via the 
hyperfine interaction. Simultaneously a nuclear spin is 
flipped. High spin flip rates can be achieved as no time 
is lost in Landau-Zener sweeps and bidirectional feedback 
is possible, thereby eliminating the feedback-free reverse 
pumping necessary for the ST + scheme. These high spin 
flip rates promise further narrowing of the nuclear field 
distribution. Our calculations indicate that dephasing 
times of up to T 2 * = 6 ps are possible with the EDSR 
scheme. 

The remainder of the paper is organized as follows. 
In section U we first briefly review the nuclear dynamics 
and then present the spin diffusion model relating pump 
rate and noise. Section El compares the results of the 
model to experimental data. EDSR and the new EDSR 
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scheme for a single quantum dot are described and the 
achievable narrowing is evaluated in section IIIII Section 
IIVI shows how the EDSR scheme can be applied to a 
double dot system with and without exchange coupling. 
A summary and discussion of results is given in section 

El 


I. NUCLEAR DYNAMICS AND SPIN 
DIFFUSION MODEL 

A. Nuclear spins 

The nuclear spins of the host material and the elec¬ 
trons in the dot interact via the Fermi contact hyperfine 
interaction described by 

Hu = I / 0 ^A Q(i )|V’i(R 1 )| 2 Se-I l (1) 

i 

as illustrated in figure [Q Here, A a (q is the microscopic 
hyperfine coupling for the different nuclear species a(i), 
which is weighted by the local electron density |^j(Ri)| 2 , 
and uq is the volume of the primitive unit cell containing 
two nuclei (in III-V materials). S e is the electron spin 
and I i the angular momentum operator of the Ith nuclear 
spin. 



FIG. 1. Illustration of the interaction between nuclear spins 
(small arrows) and the electron spin (big arrow) in the dot. 
If the wave function of the nuclei and the electron wave func¬ 
tion overlap, a finite probability of a simultaneous electron 
and nuclear spin flip occurs due to hyperfine contact interac¬ 
tion. The spatially varying polarization of the nuclei can be 
described by the function P(r). This polarization influences 
the electron’s Larmor frequency. 

In consideration of the large number of nuclear spins 
interacting significantly with the electron (typically 10 6 ), 
we describe their effect in terms of a classical effective 
magnetic field, the Overhauser field, 

B nuc = ^-^|^(R i )| 2 <I i ), (2) 

9>b ^ 

where g* ~ —0.4 is the effective g-factor, pe the 
Bohr magneton and A = 90 peV the hyperfine coupling 
strength. For simplicity, we assume a constant hyper¬ 
fine coupling strength for all species. The model can 
be extended straightforwardly to drop this assumption, 
but this would not change our overall conclusion. Due 


to the high temperature of the nuclear spin bath, the 
Overhauser field fluctuates around its average. This 
field fluctuation can be modelled as a classical, Gaus¬ 
sian distributed noise variable with a standard deviation 
of « 4mT in GaAs quantum dots, where B max is 
the maximum Overhauser field corresponding to full po¬ 
larization and N the effective number of nuclei in one 
dot— For JV 1 the Overhauser field can be expressed 
as 

B nuc = H max y d 3 r|V>(r)| 2 P(r,f). (3) 

Here P(r, t) is a dimensionless local nuclear polariza¬ 
tion density and H max = where I denotes the to¬ 

tal nuclear spin quantum number. In the following only 
the component of the Overhauser field aligned with the 
external magnetic field is evaluated because this is the 
principal component seen and influenced by the electron. 
The relative polarization density can then be expressed 
as P(r, t) = (If) /I. 

B. Spin diffusion model 

For t > 100 ps nuclear dipole-dipole interactions lead 
to a diffusion-like redistribution of the local nuclear po¬ 
larization inside and outside the quantum dotj 21 i 22 At 
magnetic fields exceeding 20 mT, hyperfine mediated dif¬ 
fusion can be neglected.— We neglect nuclear spin-lattice 
relaxation as well because the relaxation time T™ c ss 
20 min is much longer than the typical correlation time of 
the Overhauser field (seconds to about 1 minute) i 17 i 24 We 
also neglect the Knight shift, which slightly suppresses 
diffusion in the region occupied by the electron. 

The nuclear dipole-dipole interaction can thus be de¬ 
scribed by the diffusion equation^ on distance scales 
larger than the lattice spacing and time scales longer than 
the nearest-neighbour dipole-dipole interaction time. 
The corresponding dynamics of the relative polarization 
are given by 

d t P(r,f) = DV 2 P(r,f), (4) 

where D is the diffusion constant— The random nuclear 
fluctuations can be incorporated via Langevin forces. A 
spatial Fourier transformation of equation fi} gives 

d t P k (t) = -Dk 2 P k (t) + E k (t) (5) 

where is the Langevin force acting on the Fourier 
mode with wave vector k and satisfies 

( Z k (t')Z k \t ")} = S 3 , k , k .5 k , k .6(l' - t"). (6) 

Here, Ss,k,k' is the magnitude of the frequency indepen¬ 
dent noise of the Langevin forces for each Fourier mode— 
To obtain the dynamics of the total Overhauser field, 
we rewrite equation <[5|) in terms of the contributions 
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of the Fourier components to the total Overhauser field 

B nno = J2k B nuc wi th 

Bnuc (t) = fk-P k (t), ( 7 ) 

fk = 73 ma x|V’lfe and denoting the Fourier transform 
of \ip\ 2 - This transformation leads to 

k k k 

( 8 ) 

If DNP is applied, a term describing the resulting 
change in polarization has to be added. Pumping in¬ 
troduces a change of B nuc by locally polarizing spins. If 
feedback is used, the pump rate dB nuc /dt is a function of 
instantaneous B nuc . We linearise this relation by intro¬ 
ducing the feedback gain 7 pU mp = — d ^ d ^' xc / dt '> and thus 
obtain 


where Sg = Y/kfk^s,k is the spectrum of the rate of 
change of the Overhauser field. Consequently the Over¬ 
hauser field variance is directly related to this measure of 
the strength of the Langevin forces and the feedback gain 
determined by the pump rate. Note that as a result of 
neglecting the fc 2 -terms, this result is of the same form as 
previous models not explicitly considering spin diffusion, 
but with a different expression for Sg. 

The feedback gain 7 pU m P can be measured (e.g. see 
section ED or computed from a microscopic model. To 
experimentally obtain Sg, we relate it to the fluctuation 
spectrum of B nuc without DNP, Sj 3 nu( . (cu). which can be 
measured independently. Starting from equation © and 
following Ref. [26|, one obtains 


Sf 


d2 

3 (w) = max 


dt 


N J (1 + tD/aly) y/l + tD/al 


( 12 ) 


J t X B nuc(t) = -7pump X B nnc(t) 
k k 

-E^uc (t)£>fc 2 -X>S*(i). (9) 

k k 

As the effective nuclear field is weighted by the Fourier 
transformed electronic charge density \i/\ 2 , modes with a 
wave length small in comparison with the spread of the 
electron wave function are not relevant. The dominant 
low-fc modes on the other hand are slow compared to the 
feedback time scale since typical decorrelation times of 
nuclear spin diffusion are seconds , whereas feedback acts 
in the sub-second regime^ii Thus the term B k uc D k 2 

will be small compared to 7 pump and can be neglected for 
all k-modes. Finally the change in the nuclear magnetic 
field can be described as 

^B nuc (t) = -7p U mp£nuc(£) - ^2 ( 10 ) 

k 

Solving this first-order inhomogeneous linear differential 
equation in B nuc (t) leads to 

t 

Bnnc{t)= J 

— oo k 

The electronic dephasing time T% depends on the vari¬ 
ance of the Overhauser field, 

t t 

{B 2 nc (t)) = J J dt'dt”e~' yputnp( ' t ~ t '' > e~ Jpun ' p( ' t ~ t "' > 

—oo — OO 

J2fkfk'{z k (t')z k \n)- 

k,k' 

Substituting (E k (t’)E k ' (t")) from equation © and car¬ 
rying out the integrals, we thus find that the variance in 
the presence of pumping is given by 

<£nucW>pum P = 7T^-, (H) 

/pump 


for a Gaussian electron wave function |^(r)| 2 = 

1/ [(27r) 3 / 2 cr 2 a z ] exp [- (r \lo\ + z 2 / ct 2 ) /2] with a z as 
the vertical and cr± as the lateral spread of the wave 
function. Approximating the ^-dependence of the wave 
function as Gaussion is not expected to introduce ma¬ 
jor errors because the exact results for all relevant quan¬ 
tum well potentials (triangular, rectangular) also exhibit 
a quadratic maximum and a fast decay of the wave func¬ 
tion, thereby leading to the same qualitative diffusion 
behaviour. — The resulting noise correlator can be read 
off as 


(. B nuc (t)B n 


= (0))e q = %^ 


(l + tD/a 2 y) + tD/a 2 z 


( 13 ) 


On time scales shorter than the correlation time t c = 
2 

min(-^-) « Is, the correlator can be approximated to 
first order as 


{B n uc(t)B n 


:(0)) eq = %^ 


[1 - {D/a, 


2 

xy 


= ( b LM) ~ 2 Set■ 


D/2 a 2 ) t] 
( 14 ) 


The last form is obtained by again neglecting the k 2 
terms in equation ©, which is valid for large ui corre¬ 
sponding to short times t. The reason for this simple 
result is that different fc-modes have Lorentzian spectra 
with different corner frequencies l/£ c , but on the feed¬ 
back time scale (shorter than correlation time t c ), all rel¬ 
evant modes show an universal 1/w 2 behaviour leading to 
a linear time dependence of the time domain correlator. 

Note that since the different diffusion components have 
different relaxation times, we cannot relate the narrowing 
of the Overhauser distribution to the maximal achievable 
field by DNP.II In the diffusive case considered here, this 
maximum field depends on dimension, relaxation rates 
and diffusion constants. Instead, we employ the directly 
measurable fluctuation spectrum without DNP. 

While we have only considered a single quantum dot so 
far for simplicity, our model is also applicable to double 











4 


quantum dots. For the most widely studied S-To qubits, 
the nuclear field B nuc is replaced by the nuclear field 
gradient between the two dots, A B nuc = B^ uc — B^ uc . 
As one does not expect correlations between the dots, 
the total spectrum is the sum of those from each dot. 


II. COMPARISON WITH EXPERIMENT 

We now apply our model to the experiment from 
Ref. 0 performed on a gated double quantum dot in a 
GaAs/AlGaAs heterostructure. DNP was obtained from 
Landau-Zener sweeps across the crossing of the singlet 
state S and the m = 1 triplet state T + , where spin an¬ 
gular momentum can be transferred between electrons 
and nuclei. Feedback was introduced by preceding each 
sweep with a free evolution of the qubit driven by the 
current hyperfine field gradient AB„ UC . This precession 
modulates the pump rate as a function of A B* uc because 
only one of the two possible qubit states at its end allow 
spin transfer to the nuclei in the subsequent sweep. 

The independently measured noise correlator of the 
Overhauser field fluctuations in the same sample is shown 
in figure [2j By fitting the noise correlator with equation 
(fT3l) . the noise magnitude Sg can be calculated according 
to equation (ThH) yielding = 0.53mT 2 /s. u z /a xy was 
fixed at 0.32, and the fitted value of (R 2 uc) 1 ^ 2 = 2.5 mT 
agrees with typical values for T 2 *. With N = 4.3 • 10 6 
nuclear spins corresponding to <r z = 6 nm, the fit gives a 
value for the diffusion constant of D = 5 • 10~ 14 cm 2 /s, 
which is also reasonable ; 17 ’ 28 

The feedback gain 7 pU mp can be derived from the DNP- 
induced rate of change of the nuclear polarization, which 
was reported as 40mT/s in Ref. @. As the feedback 
pulse was only applied for 30 % of the time, one ob¬ 
tains a mean pump rate of Fp^ mp = 12mT/s. For the 
particular feedback implementation considered, the mod¬ 
ulation of this rate with the probability of the free evolu¬ 
tion ending in an S state leads to d ^" uc (A B* uc , Tfb) = 

r pum P ’ (cos(w L Tfb)/2 + i) with w L = 9<,MB fi Bnuc be¬ 
ing the electron Larmor frequency and tr, the dura¬ 
tion of the precession between singlet and triplet state 
which leads to the feedback effect. Hence a feedback 
gain of 7p U mp = dB nuc /dB nuc = 6.33 s^ 1 is obtained 
with Tfb = 30 ns. Equation CEB then leads to an RMS 
width of the narrowed nuclear field gradient distribution 
of (B 2 uc )pum P = 0.20 mT, corresponding to a dephasing 
time of T 2 = 197 ns. This result is a factor 2 larger than 
the measured value of T 2 = 94 ns. The discrepancy likely 
arises partly from the fact that polarization and measure¬ 
ment pulses were each applied for 30 to 40 ms, which is 
too long for the approximation of an average pump rate 
as used here to be very accurate. In addition, the actual 
pump rate during feedback may have been lower than in 
the measurement used here. 

The noise spectrum is mainly determined by dot geom¬ 
etry and independent of the feedback scheme. Therefore 



At (S) 

FIG. 2. Measured noise correlator of the Overhauser field 
fluctuations in a GaAs double quantum dot without DNP ob¬ 
tained from a time series of measurements of AB^ C via FID 
measurements with a sampling rate of 0.6 Hz. We plot the dif¬ 
ference between the correlator and the variance because this 
subtraction cancels a relatively large statistical uncertainty 
in the variance arising from the long correlation time. The 
short-time behavior of the correlator is not affected by this 
uncertainty. The slope near At = 0 yields the noise mag¬ 
nitude Sg = 0.53mT 2 /s, which can be used to predict the 
variance of the Overhauser field distribution. 


to increase the effectiveness of the narrowing of the dis¬ 
tribution, the feedback gain and accordingly the pump 
rate should be maximized. 


III. POLARIZATION SCHEME WITH EDSR ON 
A SINGLE QUANTUM DOT 

We propose a new polarization scheme, the ’EDSR 
scheme’, based on the work of Laird et al — to improve 
the feedback efficiency. Hyperfine-mediated EDSR can 
be used to control the nuclear magnetic field in the quan¬ 
tum dot as spin angular momentum conservation requires 
a simultaneous nuclear and electron spin flip. An oscil¬ 
lating electric field shifts the electron wave function in 
the dot, which in turn creates an effective perpendicu¬ 
lar ac magnetic field because the electron experiences an 
oscillating hyperfine field B^ nc 2& B„ uc drives the elec¬ 
tron spin with a Rabi frequency Dr which depends on 
the amplitude of the electric field and the hyperfine cou¬ 
pling strength. Competing flipping mechanisms such as 
spin-orbit coupling decrease the nuclear spin flip proba¬ 
bility as no nuclear spin is required to flip. However, by 
aligning the external magnetic field perpendicular to the 
direction of electron motion, the spin orbit channel can 

be suppressed^ 

To stabilize the Overhauser field, two ac electric fields 
with frequencies each slightly detuned from the desired 
electron Larmor frequency are sequentially applied to 
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FIG. 3. (a) Illustration of the pump curve created by 

the negatively detuned microwave pulse driving the transi¬ 
tion It) —► |4_) (part 1 of the EDSR feedback scheme, red) and 
the pump curve created by the positively detuned microwave 
pulse driving |4_) —> |t) (part 2 of the EDSR feedback scheme, 
blue), (b) Pump curve created by combining the positively 
and negatively detuned part of the EDSR scheme. In green 
(black) the maximum-gain (most stable) pump curve yielding 
a feedback gain 7 pump = 1.2 ■ 10 4 s _1 (7 pump = 0.8 ■ 10 4 s _1 ) 
is depicted. The parameters fulfil the conditions ^ = 1.4 

(^- = 3) and tmwSu: = 2.3 (tmw<5w = tt) and the ex¬ 
act parameters used are Hr = 1.8 ps -1 , tmw = 0.91 ps 
(tmw = 0.58ps), 5lo = 2.52 ps -1 (Six = 5.40 ps -1 ) and 
T cyc = tmw + 250 ns. The ratio of the distance between the 
pair of void (filled) diamonds to the fluctuation width is 12 
(19) so that the feedback is stable. 


drive transitions between the electron up (|t)) and the 
electron down (|4_)) state in a single quantum dot. First 
a |t) state is initialized and a microwave pulse of duration 
tmw at a frequency negatively detuned by — dui from the 
targeted electron Larmor frequency wl = g*fj ,b Bext ^ 5imc 
is applied. If the nuclear magnetic field decreases and 
reaches the value corresponding to this detuned fre¬ 
quency, the transition from |t) to |4_) is driven. Together 
with the electron spin a nuclear spin is flipped via the 
hyperfine interaction and the nuclear magnetic field in¬ 
creases (part 1). The rate of change of the nuclear field 
as a function of the nuclear field (’pump curve’) is illus¬ 
trated in figure [3^,. Next the qubit is initialized in a |4-) 


state and a microwave pulse at a frequency positively de¬ 
tuned by +<5w is applied (part 2). Now the transition 
from |4.) to |t) is driven and the nuclear field decreases 
if the field strength has reached the value corresponding 
to this positively detuned frequency. The pump curve is 
also shown in figure [3jr. Alternating between both parts 
of the scheme leads to a total pump curve with a stable 
point at the desired nuclear field wl as shown in figure [3}j. 
This dependence of the polarization rate on the nuclear 
state allows the qubit to act as a feedback loop. The neg¬ 
ative feedback reduces fluctuations by driving the system 
to the stable point. 

In order to estimate the effectiveness of the EDSR 
scheme, i.e. the narrowing of the Overhauser distribu¬ 
tion, we calculate the pump curve and thus the feed¬ 
back gain. Firstly the total spin flip probability is 
evaluated following Ref. [l9|. Because of the parabolic 
confinement potential, the displacement of the electron 
wave function is approximately proportional to the elec¬ 
tric field. Furthermore, since the displacement is much 
smaller than the spread of the electron wave function, 
the change in nuclear field is approximately proportional 
to it. We assume that for each electron spin also a 
nuclear spin is flipped. Hence the probability distribu¬ 
tion of the Rabi frequency arising from the 2D Gaussian 
distribution of the transverse nuclear field is given by 
p(D) = ^D/Dj^) exp(—D 2 /D|J, where Dr is the typi¬ 
cal Rabi frequency. Multiplying the distribution with 
the EDSR spin-flip probability from an initial |4_) state, 
we arrive at a total averaged spin flip probability for 
hyperfine-mediated EDSR of 


Pi 


A w 


tmw^r = / dx- 


2x 3 


.2 / Auj/2 \' 




Sill 


f Aw/2\' 

ww 




tmw^rJ e 


(15) 


Here Aw is the detuning of the microwave from the 
current Larmor frequency, tmw the pulse duration and 
x = the Rabi frequency relative to the typical Rabi 
frequency. We see that the pump curve depends on the 
detuning and pulse duration relative to the typical Rabi 
frequency and tmwHr, implying that the pump curve 
and in turn the narrowing is independent of the absolute 
value of Dr. The probability p-f to flip from ||) to |4.) is 
given by the same expression. Converting this spin-flip 
probability to a change in nuclear field and adding up 
the negatively (part 1) and positively (part 2) detuned 
pulses results in a pump curve of 


dBnuc ( ^nuc 


dt \ Dr 


, tmwDr = 


1 




^nuc 4” did 

Pi l --, TmwDr 


Dr 


- 


Tcyc (tmw) A 
(^nuc did 

Pi I —f;-,tmwHr 


Dr 


( 16 ) 
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Here r cyc is the duration of the feedback cycle and 
w nuc ± Sco the combined detuning of the microwave and 
the current nuclear field from the desired electron Larmor 
frequency. Maximizing the feedback gain with respect to 
Sco and tmw for Rabi frequencies in the order of MHz (as 
measured in Ref. 0 ), gives the conditions 

Sco 

—— = 1.4 and tmwScj = 2.3. (17) 

Hr 

The resulting pump curve is plotted in figure [3] in green. 
For T C y C = tmw this pump curve yields the optimal so¬ 
lution for all sets of parameters Sui, tmw and Hr ful¬ 
filling (1171) but if r cyc > tmw due to the initialization 
time, the effective feedback gain still depends slightly on 
the individual parameters and is not completely fixed by 
the conditions. Using the largest value from Ref. |l9], 
Hr = 1.8 ns -1 , for the Rabi frequency leads to the pa¬ 
rameters tmw = 0.91 ps and Sco = 2.52 ps -1 and a feed¬ 
back gain of 7 pum p = 3.5- 10 4 s -1 . T cyc is typically 250 ns 
longer than tmw due to the initialization process. 

In combination with the noise magnitude from figure 
[2j the above feedback gain leads to an RMS width of 
the narrowed nuclear field distribution of {B 2 ^) 1 / 2 = 
6.6 • 10 -3 mT (equation 1111) ) and a dephasing time of 
T 2 * = 6 ps. Thus, the proposed pump scheme promises an 
improvement of dephasing time by more than two orders 
of magnitude in comparison to the intrinsic T 2 * = 15 ns. 
The nuclear field gradient variation in turn is larger by a 
factor of \[2. 

For the feedback to be stable, i.e. to prevent the mean 
of the distribution from escaping, the nuclear field fluctu¬ 
ations have to be at least one order of magnitude smaller 
than the the distance of the global maximum and mini¬ 
mum of the pump curve (see distance between ♦ and 0 in 
figure [3]). Although the narrowing is independent of Hr, 
a smaller Hr deteriorates the feedback stability. Further¬ 
more, it may be advantageous if the second order local 
maxima and minima are small, i.e. the pump curve has 
to be smooth. Appropriate parameters for the latter are 
= 3 and tmw& = 7r, which give a dephasing time 
of T 2 = 5ps with Hr = 1.8 ps -1 , tmw = 0.58 ps and 
Sco = 5.4 ps -1 . The respective curve is shown in figure [3] 
by the black trace. The ratio of peak distance to fluctu¬ 
ation width is 12 for the maximum (green) curve while 
it is 19 for the smoother (black) pump curve. Thus the 
stabilization criterion is fulfilled for both curves though 
the smoother curve is more stable. 

So far we have only considered polarization noise aris¬ 
ing from the diffusive spin dynamics. In addition, shot 
noise generated by the discrete flipping of spins is a noise 
source intrinsic to the feedback scheme. The spectral 
density of this shot noise to be added to the noise from 
diffusion is given by 


flip. The remaining factor reflects the variance of the bi- 
nomially distributed number of flips. The flip probability 
at the locking point is pm P = 0.41 for the curve leading 
to the largest pump rate. With an average cycle time 
of T cyc = 1 ps, -B max = 4 T and N = 4 • 10 6 the noise 
magnitude is 5 sn = 0.24mT 2 /s. This value is compara¬ 
ble to the diffusion noise. 5 sn has to be added to Sg in 
equation m • Consequently, T 2 is reduced by a factor of 
order unity. This shot noise can be reduced by a factor 
of order unity by choosing a smaller Hr corresponding 
to a larger T cyc and a reduced microwave power with¬ 
out changing the feedback gain 7 pU mp- Ultimately, this 
reduction is limited by the distance between the peaks 
of the pump curve becoming comparable to the fluctu¬ 
ation width, which will compromise the stability of the 
feedback scheme. 


IV. GENERALIZATION TO A DOUBLE 
QUANTUM DOT AND INFLUENCE OF 
EXCHANGE COUPLING 

A. Decoupled electrons 

The EDSR pump scheme described for a single dot in 
section mu can easily be adapted for a double dot system 
for which the spin states are |LR), where L, R= {'}, |} 
refer to the electron state in the left and right dot, re¬ 
spectively. In this case, the individual magnetic fields in 
the dots are stabilized so that the field gradient AB^ UC is 
stabilized as well. Initialization of the qubit in |4_t) (I'M-)) 
can be achieved by adiabatic preparation from a singlet 
5(0,2), where both electrons are in the same dot, (and 
applying a 7r-pulse). Consequently both dots are jointly 
initialized. Then the microwave pulses from section ED 
are applied to each individual dot which can be addressed 
either locally or through frequency selection. For the lat¬ 
ter the desired nuclear field gradient A B* uc has to be 
large enough so that |j_t) and |t4-) are energetically well 
separated. For the right dot the transitions driven by 
the microwave are illustrated in figure [4p,. The electron 
is driven into the triplet state T + = |4_|) (T_ = |tt)) with 
a positively (negatively) detuned microwave pulse if the 
nuclear field in the right dot has fluctuated downwards 
(upwards). Consequently the nuclear field in the right 
dot increases (decreases) dependent on the state of the 
nuclear field in that dot. The same scheme is applied to 
the left dot at a different frequency, either simultaneously 
or sequentially. A sketch of the pulsing scheme is shown 
in figured}:). 


B. Exchange coupling 


5sn = 


Pfi ip * (1 Pflip) (B : 


1 eye 


N 


(18) 


where is the magnetic field change per nuclear spin 


In a double dot system, an alternative possibility to 
drive transitions is the exchange coupling J{e)^ At fi¬ 
nite J the microwave pulse leads to an oscillation of 
the exchange coupling described by J(t) ps J(eo) + 
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requiring J(e) < g* e \i bAB^ uc . On the other hand, the ex¬ 
change energy has to be large enough to drive the tran¬ 
sitions as the level splitting becomes J-independent for 
J -C g*/isAi?^ uc and SJ is bounded by J(e). Insert¬ 
ing typical experimental parameters for S — To qubits in 
GaAs (SJ/h = 50 MHz, AB^ UC = 5mT, Bf xt = 500 mT) 
yields a Rabi frequency of Hr = 0.81 ps _1 . The Rabi fre¬ 
quency is of the same order of magnitude as the EDSR 
Rabi frequency. However, a spin-orbit contribution to 
the S — T + coupling with a direct e-dependence^ may 
introduce another strong driving channel in this regime, 
which would have to be avoided carefully. 


FIG. 4. (a) Energies of the relevant double dot spin states 
as a function of the detuning in the (1,1) charge state for 
A-Bn UC > d(e). The detuned microwave signals are indicated 
by arrows. A drift of AB^uc brings one closer to and the other 
further from resonance, (b) Schematic of e(t) for the posi¬ 
tively detuned (top) and negatively detuned (bottom) feed¬ 
back pulses. 


^j-Se cos(wmw 0 = >/(eo) + cos(wmw^) where Se is the 
amplitude of the detuning oscillation dependent on the 
power of the EDSR burst, which in turn gives the ex¬ 
change oscillation amplitude 5J. This exchange oscilla¬ 
tion can be used as an alternative nuclear spin flipping 
mechanism. 

Without EDSR the transitions | S) —> \T±) are driven 
by the perpendicular magnetic field difference AB^ UC , but 
as ABi c < H| xt the transition probability is negligible 
due to conservation of energy. This suppression is over¬ 
come by the ac-drive, which corresponds to a Hamilto¬ 
nian with a time-dependent variation of the level splitting 
given by 


H = 


hLO z (t) 

-^- a z + 


Huj± 

2 ay ’ 


(19) 


where cr, are the Pauli matrices, huj z = — J(t) 

and hoj _l = ^ggB^nuc* Transforming the eigenbasis for 
static e (Se = SJ = 0) and neglecting the time depen¬ 
dence of the level splitting in this basis leads to a stan¬ 
dard Rabi Hamiltonian with a Rabi frequency of 


n 


J _ 
R — 


+ w 2 


SJ_ 

4ft' 


( 20 ) 


For a typical external magnetic field R| xt AB I ) UC , the 
Rabi frequency can be approximated to 


^R = 


AB 


_L 

nuc 


B? xt - J (e 0 )/(g*v B ) 


5J_ 
4 ft' 


( 21 ) 


In order for the exchange flipping mechanism to cre¬ 
ate a feedback effect, the two eigenstates need to have 
substantially different components of |t4-) and |4_t), thus 


V. SUMMARY AND CONCLUSION 

In conclusion, we have developed a model that pro¬ 
vides a simple relation between the feedback gain of DNP 
procedures with feedback, the spectrum of nuclear spin 
fluctuations, and the achievable degree of narrowing of 
the Overhauser field distribution. The dynamics of the 
nuclear spins are modelled as a diffusive process. The 
results are in agreement with earlier experiments. We 
then proposed a new feedback scheme for gated quan¬ 
tum dots based on EDSR which should achieve signifi¬ 
cantly larger pump rates and thus feedback gains. Our 
model predicts dephasing times of up to T 2 * = 6 ps for 
this scheme when considering only the diffusive fluctua¬ 
tions of the z-component of the Overhauser field. Inter¬ 
estingly, this value is comparable to the measured value 
of the Hahn-echo coherence time, which was found to be 
on the order of 30 ps at sufficiently high magnetic fields£L 
Achieving such long dephasing times would significantly 
enhance the qubit control fidelity. In practice, one may 
however find that effects not considered here impose fur¬ 
ther limitations. In particular, the faster dynamics of 
transverse nuclear field components have pronounced ef¬ 
fects on Hahn echo amplitudes at applied fields of a few 
100 mT and below, and are expected to affect FID mea¬ 
surements even more^i Considerably shorter intrinsic 
FID times have been predicted when taking transverse 
terms into account. 21 - To partially overcome this limita¬ 
tion, it may turn out important to synchronize both the 
feedback scheme and measurements to the relative Lar- 
mor precession of different spin species. On the other 
hand, it is conceivable that even better results can be 
obtained at high fields because of an expected cut-off 
of the spin diffusion spectrum at frequencies exceeding 
the dipolar nuclear coupling strength. In any case, ex¬ 
perimentally probing the limits of the proposed EDSR 
feedback scheme is likely to yield new insights into elec¬ 
tron spin dephasing by the nuclear spin bath as it may 
give access to the so called intrinsic FID, which is so far 
a purely theoretical construct ; 21 ! 32 

We note T 2 * times for gated dots approaching 3 ps 
have already been demonstrated by a measurement based 
narrowing procedure^ Compared to this approach, our 
method would have the advantage of achieving a pre- 























determined center of the Overhauser field distribution 
rather than narrowing around a random value. 
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